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The mass spectrum of the 3d 3-state Potts model is considered in the broken phase (a) near the 
second order Ising critical point in the temperature-magnetic field plane and (b) near the weakly 
first order transition point at zero magnetic field. In the case (a), the mass spectrum is compared 
with the prediction from universality of mass ratios in the id Ising class; in the case (b) a mass 
ratio is determined to be compared with the corresponding one in the spectrum of screening 
masses of the {1>+Y)d SU(3) pure gauge theory at finite temperature in the deconfined phase near 
the transition. 
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1. Introduction 

The Svetitsky-Yaffe conjecture [|l|] establishes that l)-dimensional SU(N) pure gauge the- 
ories at finite temperature, which undergo a confinement/deconfinement phase transition associated 
with the breaking of the center of the gauge group Z(N) [Q, ||], are in the same universaUty class of 
J-dimensional Z(N) spin models, which undergo an order/disorder transition. This implies that the 
gauge theory and the spin model present the same universal quantities (critical indices, amplitude 
ratios, ...) at criticality. This has been verified in the case of a second order transition by numerical 
investigation in several papers (see, for instance, Refs. and, for a review, Ref. [§]). Moreover 
recently it has been conjectured that also the ratios of massive excitations in the broken phase are 
universal [§] and a numerical evidence of this conjecture has been given in the case of 3d Ising 
universality class [^, ||, 

In this work we consider the 3d 3-state Potts model, which has an interesting phase diagram in 
the inverse temperature-magnetic field plane (see Fig. |l]) and provides a test-field for universality 
ideas in two different regions: (a) near the critical endpoint in the (j3 — h) plane, belonging to the 3d 
Ising universality class [^] and (b) near the h = weakly first order transition in the broken phase. 
In both cases we calculate the ratio between two massive excitations to be compared in the case 
(a) with the corresponding mass ratio of the 3d Ising class; in the case (b) with the corresponding 
ratio of screening masses in the (3+l)d SU(3) pure gauge theory at finite temperature near the 



deconfinement transition in the broken phase, calculated in Ref. ||lip. The reason of the latter 
comparison is to check if and to what extent Svetitsky-Yaffe conjecture holds in the case of weakly 
first order transition. 



2. The 3d 3-state Potts model and massive excitations 



The 3d 3-states Potts model [|I2[ [T^ is a spin theory in which the fundamental degree of 
h 




Figure 1: Qualitative phase diagram of the 3d 3-states Potts model in the (j3,/;)-plane: the solid line in 
bold is the line of first order phase transitions. j3( is the order/disorder inverse transition temperature at zero 



magnetic field; ^ and T are the symmetry-breaking and the temperature parameters of the Ising theory (2.5) 
and (Tc^c) is the endpoint in the 3d Ising class. 
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freedom, 5,, defined in the site / of a 3-dimensionaI lattice, is an element of the Z(3) group, i.e. 

Si = e'l'"'% a,- = {0,1,2} . (2.1) 

The Hamiltonian of the model is 

^ = - 1 {s]sj + s]s^) = -/3 £ 5,,,,^. , (2.2) 

^ iU) (ij) 

up to an irrelevant constant. Here, j8 is the coupling in units of the temperature and the sum is 
done over all the nearest-neighbor pairs of a cubic lattice with linear size L. It is well known that 
this system undergoes a weakly first order phase transition [[T^, associated with the spontaneous 
breaking of the Z(3) symmetry. The order parameter of this transition is the magnetization, 

W = (^I^/)- (2.3) 

In presence of an external magnetic field it is convenient to work with an Hamiltonian written in 
terms of the a, degrees of freedom, instead of the Si ones. For a uniform magnetic field along the 
direction Ou with strength h in units of the temperature, the Hamiltonian is 

// = -/3 £ 5o,^o^ - /j£ 5o..c„ = -PE - hM , (2.4) 

where E is the internal energy and M is the magnetization. The magnetic field breaks explicitly the 
Z(3) symmetry. However, first order transitions still occur for values of the magnetic field strength 
h below a critical value he, the transition coupling decreasing with increasing h. The line of first 
order phase transitions ends in a second order critical point Pc = {pc,hc) (see Fig. [T]), belonging to 



the 3d Ising class [10]. The Hamiltonian in the scaling region of the critical point can be written 
as 

H = -rE-^M , (2.5) 

where E and M are the Ising-like energy and magnetization operators and T and ^ the corresponding 
temperature-like and symmetry-breaking-like parameters. This means that (M) is the new order 
parameter. Close enough to Pc, the following relations hold, 

M = M + sE, E = E + rM, (2.6) 

where the mixing parameters (r, s) have been determined numerically for several lattice sizes L in 
Ref. [10]. The T-direction identifies the first order line (see Fig. |l|). 



Among the quantities relevant in the description of a phase transition there is the correlation 
function of local order parameter: in our case of the 3d 3-states Potts model this is just the local 
spin Si. The point-point correlation function is defined as 



rt,{r) = {sjsi,)-{s]){si,), (2.7) 

where / and /q are the indices of sites and r is the distance between them. The large-r behavior of 
the point-point correlation function is determined by the connected correlation length of the theory. 
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^0, or, equivalently, by its inverse, the fundamental mass. It is convenient to study the connected 
wall-wall correlator in the jc-direction defined as 

= jY,{w{xo+x)^w{xo)) - (w(x + xo)"^)(w(xo)) , w(x) = jjL'^iw} (2-8) 
^ -vo ^ y,z 

where w{x) represents the spin average over the "wall" at the coordinate x. 
The general behavior for the function G{x) is 

GW=£a„e-"'"\ (2.9) 

n 

where niQ is the fundamental mass, while m\, mj, ... are higher masses with the same angular 
momentum and parity quantum numbers of the fundamental mass. Mass excitations in channels 
different from 0+ can be determined by a suitable redefinition of the wall average. The fundamental 
mass in a definite channel can be extracted from wall-wall correlators by looking for a plateau of the 
effective mass, mgff(x) = — \nG{x)/G{x— 1), at large distances. In the present work we consider 
only the 0+ and the 2+ channels; the local variables to be wall-averaged as in have been 
defined in the following way: 

^\x,y,z} (") ~ '^{-t ^^{x,y+n,z] + ^{x,y,z+n} ) i ^\x,y,z?s (") ~ ^{x,y,z} {^{x,y+n,z] ~ ^{x,y,z+n} ) ■ 

According to the conjecture of universal spectrum, ratios between massive excitations in the 
broken phase must be equal in theories belonging to the same universality class. The aim of this 
work is, then, to compare m2+/mQ+ ratio in two cases: (a) near the 3d Ising critical endpoint to 
verify whether this is consistent with the class prediction; (b) in the broken phase near the first 
order transition in absence of the magnetic field. In the first case, the procedure to determine 
the fundamental masses in the two channels of interest is the same outlined above, with the only 
difference that we need to use the correct local order parameter to build correlators. We have 
defined this local variable, m,, in such a way that it reconstructs the global magnetization operator 
M after summation over the whole lattice: 

mi = 8a,oi, + :r £ . (2.10) 

A 

Indeed, it is easy to see that M = Y^i^i. In the case (b), results are to be compared with the would- 
be universal partner of the Potts model, which is the (3-1-1) J SU(3) pure gauge theory at finite 
temperature, to test if Svetitsky-Yaffe conjecture holds. Generally universality applies in case of 
critical transitions, i.e. where correlation length diverges. In the case of /j = Potts model, as well 
as for SU(3) pure gauge theory, transition is weakly first order (<§ keeping finite); this means that E, 
takes values much larger than the lattice spacing, although finite. 



3. Numerical results 



We have performed numerical Monte Carlo simulations of the 2>d 3-states Potts model using 



a cluster algorithm [15, 16] to reduce the autocorrelation effects. In order to minimize the finite 
volume effects, we have imposed periodic boundary conditions. Data analysis has been done by 



4 



Universality and massive excitations in 3d 3-state Potts model 



M. Gravina 




Figure 2: Effective mass in the 0+ and 2+ channels as a function of the separation between walls on the 
{y,z) plane at ^ = and T=0. 37248, determined from the configurations belonging to the "right-peak" in 
the thermal equilibrium ensemble. 



the jackknife method apphed to bins of different lengths. In both the regions studied we have 
seen tunneling between degenerate minima near the transition point. This finite volume effect can 
spoil mass measurements in the scaling region and must be treated carefully. Depending on the 
order of the transition, tunneling effects show up differently and must be removed accordingly. 
In the region (a) we have performed simulations on 70^ lattices for which the mixing parameters 
appearing in (U) turn to be s{L = 70) = -0.689(8) and r(L = 70) = 0.690(3) In this case 
tunneling occurs between two Z(2)-broken minima; this effect is removed simulating the system far 
enough = t,c,T = 0.37248) from the transition, in the scaling region, where two corresponding 
Z(2)-peaks ("right" and "left") are well separated. As an example, in Fig. |^ it is displayed the 
behavior of the effective mass in the 0+ and 2+ channels for the data in the "right" Z(2)-peak; the 
fundamental masses in each channel turn to be the plateau values of the effective mass. We have 
found 



"right-peak" (stat. 115K): amQ+ 
"left-peak" (stat. 85K): amQ+ 



0.0725(63) , am2+ 
: 0.0714(40) , £?m2+ 



: 0.1981(87), ^=2.73(36); 



0.1959(80) , 



2.74(27) . 



We can see that the mass ratios are consistent, as expected and, moreover, they are compatible with 
the value of the M Ising class [^, m2+ /mQ+ = 2.59(4). 

In the region (b) we have performed simulations on 48^ lattices for several values of the cou- 
pling j8 in the broken phase of the 3(i 3-state Potts model at zero magnetic field. 

Close enough to j3f(L = 48)=0.550538 [14], the scatter plot of the complex order parameter 



{S) shows the coexistence of the symmetric phase (points around (0,0) in the \m{S) - Re (5) plane 
in Fig. H) and of the broken phase (points around the three roots of the identity in Fig. ^. 

To remove the former kind of tunneling we just moved away from the region across transition 
where symmetric phase appears; this happens up to j3 = 0.5508. For the latter one, we applied an 
unambiguous rotation of all configurations belonging to the complex sectors to that along the real 
axis; this allows to improve statistics. We performed calculation of fundamental masses in 0^ and 
2+ channels, as done in the case (a), up to j3 = 0.6. In Fig. ^ we have plotted the behavior of mo+ 



and m.2+ versus jS (for more details look at [18]). We have determined the ratio m2+ /niQ+ for several 
j8 values in the region [j8f,0.56] which is plotted in Fig. ^. This ratio remains practically constant 
in the considered region, this suggesting that the correlation length (§2 associated to the channel 2+ 
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Figure 3: (Left) Scatter plot of the complex order parameter S at the transition point j3,=0. 550538 at 
zero magnetic field. Both the symmetric and the broken phases are present. (Right) Fundamental masses in 
the 0+ and 2+ channels as functions of j3, in the broken phase near j3( (vertical line). 



(i§2 = l/m2+) scales in the same way of the fundamental one (i^o = l/'^o+)- We can take as our 
estimation of the mass ratio the value 



m2+ 



2.43(10) 



(3.1) 



determined by taking value and error of the point with the smallest error belonging to the plateau. 
In Fig. ^ we plotted the corresponding ratio in SU(3) pure gauge theory [|ll|] in the broken phase 
near the deconfinement transition and compared them with the result for the Potts model. As we 
can see the screening ratio turns to be larger than the corresponding one in the spin model. 




.66 5.68 5.7 5.72 5.74 S.76 5.78 5.8 5.82 5.84 5.86 5.88 5.9 5.92 5.94 



Figure 4: (Left) m2+ ()3 ) /mQ+ (/3 ) for (5 varying in the scaling region. (Right) Ratio m2+ / mQ+ as a function 
of the coupling j3 in the {3+l)d SU(3) pure gauge theory near the deconfinement transition in the broken 
phase [pi]]; the upper horizontal line is the fit of masses closest to transition and consistent with a constant; 
the lower one is the fit of the corresponding mass ratio in 3d 3-state Potts model. 



4. Conclusions 

In this work we have studied massive excitations of the 3d 3-states Potts model near the Ising 
critical point on the inverse temperature - magnetic field phase diagram and near the transition point 
at zero magnetic field. 



6 



Universality and massive excitations in 3d 3-state Potts model 



M. Gravina 



We have found evidence that the mass ratio m2+ /mo+ near the Ising critical point is compatible 
with the prediction from universality, thus supporting the conjecture of universal spectrum. 

In the broken phase of the scaling region near the transition in absence of the external source, 
we have found m2+ /m()+=2.43(10). This ratio turns to be ~30% larger than the ratio of the lowest 
massive excitations in the same channels of the (3+l)J SU(3) pure gauge theory at finite tempera- 
ture in the broken phase. This can be taken as an estimate of the level of approximation by which 
the Svetitsky-Yaffe conjecture, valid in strict sense only for continuous phase transitions, can play 
some role also in this case of weakly first order transition. 
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